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Approximating  Conductive  Ellipsoid  Inductive  Responses 
Using  Static  Quadrupole  Moments 


Torquil  Smith 

Lawrence  Berkeley  National  Laboratory 
Berkeley,  California  94720 

Introduction 

Smith  and  Morrison  (2006)  developed  an  approximation  for  the  inductive  response 
of  conducting  magnetic  (permeable)  spheroids  (e.g.,  steel  spheroids)  based  on  the  induc¬ 
tive  response  of  conducting  magnetic  spheres  of  related  dimensions.  Spheroids  are  axi¬ 
ally  symmetric  objects  with  elliptical  cross-sections  along  the  axis  of  symmetry  and  cir¬ 
cular  cross  sections  perpendicular  to  the  axis  of  symmetry.  Spheroids  are  useful  as  an 
approximation  to  the  shapes  of  unexploded  ordnance  (UXO)  for  approximating  their 
responses.  Ellipsoids  are  more  general  objects  with  three  orthogonal  principal  axes,  with 
elliptical  cross  sections  along  planes  normal  to  the  axes.  Ellipsoids  reduce  to  spheroids  in 
the  limiting  case  of  ellipsoids  with  cross-sections  that  are  in  fact  circles  along  planes  nor¬ 
mal  to  one  axis.  Parametrizing  the  inductive  response  of  unknown  objects  in  terms  of  the 
response  of  an  ellipsoid  is  useful  as  it  allows  fitting  responses  of  objects  with  no  axis  of 
symmetry,  in  addition  to  fitting  the  responses  of  axially  symmetric  objects.  It  is  thus 
more  appropriate  for  fitting  the  responses  of  metal  scrap  to  be  distinguished  electromag- 
netically  from  unexploded  ordnance.  Here  the  method  of  Smith  and  Morrison  (2006)  is 
generalized  to  the  case  of  conductive  magnetic  ellipsoids,  and  a  simplified  form  used  to 
parametrize  the  inductive  response  of  isolated  objects.  The  simplified  form  is  developed 
for  the  case  of  non-uniform  source  fields,  for  the  first  eight  terms  in  an  ellipsoidal  har¬ 
monic  decomposition  of  the  source  fields,  allowing  limited  corrections  for  source  field 
geometry  beyond  the  common  assumption  of  uniform  source  fields. 
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Ellipsoidal  Coordinates 

The  equation  for  an  ellipsoid  with  semi-axes  a’,  b’,  and  c’  in  the  x,  y  and  z  direc¬ 
tions  is 


2  2  2 
y  z  . 

— H - —  1  • 
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'1 
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^'2 


The  equation  for  ellipsoidal  coordinates  ^2^  and  <^3  is 


2  22 
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with  b  <  a  hy  convention  (Morse  and  Feshbach,  1957,  or  Hobson,  1931,  with  roles  of  x 
and  z  interchanged).  Equation  (2)  is  a  cubic  equation  in  when  v,  y,  and  z  are  held 
fixed,  so  has  three  solutions  The  three  coordinates  <^i,  ^2^  and  1^3  are  such  that 
-b  <  ^2,  <  b  <  1(^2^  <  a  <  \^i\.  Solutions  in  x,  y,  and  z  of  Equation  (2),  for  >  a  held  fixed, 
define  surfaces  of  constant  which  are  confocal  ellipsoids  with  short  axis  in  the  x  direc¬ 
tion  and  long  axis  in  the  z  direction.  That  is,  they  give  solutions  of  Equation  (1)  with 

c'  =  a  =  -  a^,  b'  =  -  b^,  where  a  <  b'  <  c' .  For  computations  involving  a 

conducting  ellipsoid  with  semi-axes  a,  b',  and  c,  the  ellipsoidal  coordinate  =  c  delim¬ 
its  the  conducting  ellipsoid,  with  values  a  <  \^y\  <  c  defining  interior  ellipsoidal  sur¬ 
faces  and  ^y  values  c'  <\^y\  defining  exterior  ellipsoidal  surfaces.  The  limiting  case  ellip¬ 
soid  with  ^y  =  ±a  corresponds  to  a  flattened  oval  in  the  y-z  plane  with  zero  thickness  in 
the  X  direction.  Holding  one  of  the  other  two  ellipsoidal  coordinates  ^2  or  ^3  constant 
defines  (hyperboloid)  surfaces  which  are  orthogonal  to  each  other  and  to  the  ellipsoids  of 
constant  ^y .  Coordinate  ^2  has  two  branches  joined  at  ^2  =  and  ^2  two  joined  at  1^3  = 
(Hobson,  1931).  An  ellipsoidal  coordinate  system  is  determined  by  the  choice  of  princi¬ 
pal  ellipsoid  axes  x,  y,  and  z,  and  parameters  a  and  b.  Coordinates  x,  y,  and  z  are  given 
by 

_  s^ify  - 

a(a^  -  b^y^^ 
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b(a^  — 


e2^1/2 


(3) 


_4|2|3 
^  ab  ’ 

where  5^  =  ±1  =  sign(x)  according  to  the  branch  of  ^2^  5^  =  ±1  =  sign(y)  according  to  the 
branch  of  1^3,  and  sign(^i)  =  sign(z). 

Ellipsoidal  coordinates  are  useful  in  electromagnetic  induction  problems  in  ellip¬ 
soidal  objects  at  frequencies  (or  times)  at  which  currents  in  the  medium  containing  the 
ellipsoid  object  can  be  neglected.  At  these  frequencies  (times)  the  magnetic  field  outside 
the  object  is  curl  free  (as  well  as  divergence  free)  so  can  be  parametrized  as  the  gradient 
of  a  potential  satisfying  Laplace’s  equation.  In  ellipsoidal  coordinates,  solutions  of 
Laplace’s  equation  are  separable;  they  can  be  written  in  terms  of  products  of  functions 
and  of  the  three  ellipsoidal  coordinates:  Ej^(^2)  or 

E^(^2)  where  Ef^(^)  and  are  Lame  functions  of  the  first  and  second 

kind,  and  m  and  p  are  indices  which  distinguish  the  different  Lame  functions  (e.g.,  Morse 
and  Leshbach,  1953).  Lame  functions  of  the  first  kind  E^{^)  are  smooth  and  bounded  for 
^  ^  a  (e.g.,  at  the  origin  of  v,  y  and  z)  and,  except  Eq{^)  which  is  constant,  increase 
without  bound  for  ^  00,  so  they  correspond  to  external  sources,  and  Lame  functions  of 
the  second  kind  Ej^(^)  increase  without  bound  for  ^  a  ,  and  vanish  as  ^  00,  so  those 
correspond  to  responses  of  the  ellipsoids. 


Sphere  Ellipsoid  Response  Approximation 

Lor  axially  symmetric  objects  (e.g.,  a  spheroids),  axial  and  transverse  (equivalent 
dipole)  polarizability  responses  m„^(t)  and  give  the  strength,  as  a  function  of  time 
(t),  of  magnetic  dipole  moments  in  axial  and  transverse  directions,  induced  by  uniform 
(primary)  magnetic  fields  of  unit  nominal  amplitude  in  those  directions,  with  a  specified 
time  function  of  primary  magnetic  field  variation  (e.g..  Smith  and  Morrison,  2004). 
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Strictly  speaking,  polarizabilities  and  should  have  dimensions  of  Amp-m^/ 

(Amp/m)  =  m  .  However,  since  induction  coils  are  sensitive  to  the  rate  of  change  of  mag¬ 
netic  field,  we  let  niaxit)  and  nijit)  represent  the  time  rate  of  change  of  polarizabilities, 
which  have  dimensions  m  /s,  but  will  refer  to  them  as  ‘polarizabilities’.  Smith  and  Mor¬ 
rison  (2006)  approximate  the  axial  and  transverse  equivalent  dipole  polarizability 
responses  of  conductive  magnetic  (permeable)  spheroids  using  the  polarizability 
responses  of  spheres  of  related  sizes; 


^ax(0  ^  ax^  sphere(c 


(4) 


^t^sphereibgg)(.^^  ■ 

Letting  c'  be  the  axial  semi-diameter  of  the  spheroid,  and  b'  be  the  transverse  semi-diam¬ 
eter  of  the  spheroid,  the  transverse  polarizability  response  of  the  spheroid  is  approxi¬ 
mated  as  proportional  to  the  polarizability  response  of  a  sphere  of  the  trans¬ 

verse  semi-diameter  (bgjf  =  b'),  and  the  axial  response  is  approximated  as  proportional  to 
the  response  m^pf^ggg^g^pt)  of  a  sphere  of  radius  Cg^,  where  for  prolate  spheroids  {c  >  b'), 

Cgff  =  c  and  for  oblate  spheroids  {c  <  b'),  Cgjf  =  ^(c'^  +  b'c')ll,  with  the  latter  formula 
found  empirically.  This  approximation  will  be  denoted  a  sphere- spheroid  approximation. 
Here  we  have  written  the  approximation  in  terms  of  semi-diameters  (radii)  instead  of 
diameters  to  be  consistent  with  semi-axes  used  to  characterize  ellipsoids  and  ellipsoidal 
coordinates.  Proportionality  constants  and  Vt  are  ratios  of  differences  of  polarizabili¬ 
ties  at  high  and  low  frequency  limits; 


^(spheroid)^  ^  _  ^(spheroid)^ 


(0=0 


^sphere(Cgff'^(o=o 
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sphere{cp  ®=0 


(5) 


L  = 


(spheroid), 

vtif  L 
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sphere(bpco=i 


-  m 


-  m 


(spheroid), 

t 


sphere(bp  ®=0 


where  the  high  and  low  frequency  limit  polarizabilities  are  given  by  Smith  and  Morrison 
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(2006)  and  are  readily  calculated.  Generalization  to  ellipsoids  is  straight  forward.  Pro¬ 
late  spheroids  are  the  limiting  case  of  ellipsoids  with  minor  semi-axis  length  approaching 
their  intermediate  semi-axis  length  {a  b')  with  z  as  their  axis  of  symmetry.  Oblate 
spheroids  are  the  limiting  case  of  ellipsoids  with  major  semi-axis  approaching  their  inter¬ 
mediate  semi-axis  length  (c'  b')  with  x  as  their  axis  of  symmetry.  We  approximate  the 

polarizability  response  of  conducting  magnetic  ellipsoid  to  uniform  magnetic  source 
fields  (of  a  given  time  dependence)  in  the  ellipsoid  principal  directions  as 

~  ^ sphere(a^ff)0'^  j 


nry(t)  ^yf^sphere{b^ff)(.^^ 


(6) 


^z(0  ^z^sphere(Cgg)(.^^  ’ 


where 


^(ellipsoid)^  ^  _  ^(ellipsoid)^ 


o)=0 


^sphere(o^ff'}m=o 


m 


sphere(Ogff)a=0 


(7) 


where  w  stands  for  one  of  x,  y,  or  z,  and  o^jf  stands  for  the  corresponding  choice  among 
Ggff,  bgff,  and  c^jf.  For  ellipsoids,  the  appropriate  effective  sphere  dimensions  are  given 
by 


/  xl/2 

a  +  a  b 


’eff- 


b\ 


'ejf- 


(8) 


where,  as  before,  a,  b',  and  c'  are  the  semi-axes  of  the  ellipsoid.  Effective  sphere  semi¬ 
diameters  (8)  have  been  chosen  to  be  consistent  with  those  used  for  prolate  and  oblate 
spheroids,  after  allowing  for  the  difference  in  symmetry  axis  between  prolate  and  oblate 
spheroid  treated  as  limiting  cases  of  ellipsoids  with  a  <b'  <  c.  Formula  (7)  for  Vy 
and  Vj.  requires  the  zero  frequency  and  high  frequency  limit  uniform  source  field 
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polarizabilities  for  the  ellipsoid  in  source  fields  oriented  in  the  x,  y,  and  z  directions 
respectively.  These  are  derived  from  solutions  of  Laplace’s  equation  in  ellipsoidal  coor¬ 
dinates  for  an  external  field  that  is  the  gradient  of  a  potential  varying  in  proportion  to  x, 
y,  or  z. 


Ellipsoidal  Harmonics  and  the  Solution  of  Laplace’s  Equation  in  Ellipsoidal  Coordinates 

As  noted  above,  solutions  of  Laplace’s  equation  in  ellipsoidal  coordinates  can  be 
expressed  in  terms  of  products  of  Lame  functions  of  ellipsoidal  coordinates, 

(9) 


for  external  sources,  and 


(10) 

for  internal  sources.  Eorms  (9)  and  (10)  are  known  as  external  and  internal  ellipsoidal 
harmonics.  The  zero  frequency  limit  induction  problem  is  solved  using  a  term  of  each 
form  outside  the  ellipsoid,  and  a  single  term  of  form  (9)  inside  the  ellipsoid,  that  is, 

d=[ao^)^(^i)  +  /lo^^(^i)]^)^(^2)£)^(^3)  ,  outside,  ,  (H) 


Letting  be  the  coordinate  of  the  ellipsoid/exterior  interface,  the  zero  frequency 

limit  problem  is  to  match  H  x  n  and  B  ■  n  at  =  i^o  »  where  n  is  the  surface  normal  vector 
(in  the  direction).  Of  course  B  =  //H  =  //Vp.  Letting  magnetic  permeabilities  be 
//  =  //q  outside  and  //  =  //i  inside  the  ellipsoid,  this  yields 


^0^(o=0  ^olffl=0  (Mr  1) 


U'(^o) 


-Mr 


E(^o)J 


(12) 


where  '  denotes  differentiation,  =  MiIMo^  and  indices  p  and  m  have  been  elided. 
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The  high  frequency  limit  induction  problem  is  the  same  as  the  low  frequency  limit 
problem  with  B  ■  n  =  0  set  at  =  i^o-  Solution  of  this  case  results  in  aoU=oo  /?oU=oo 
obeying  Equation  (12)  in  the  limit  of  0,  that  is, 


/^(^Q)=00  ^0^O)=O 


Em 

Em 


(13) 


Equations  (12)  and  (13)  solve  low  frequency  and  high  frequency  limit  problems  for  exter¬ 
nal  sources  arising  from  a  single  external  ellipsoidal  harmonic  term. 

Eame  functions  E^{^)  and  E^{^)  depend  on  coordinate  system  parameters  a  and  b, 
which  are  constant  for  a  given  ellipsoidal  coordinate  system,  so  are  omitted  as  arguments 
of  and  F^.  Eollowing  Morse  and  Eeshbach  (1957),  the  first  few  are 


2^1/2 


2\V2 


,(14) 


'2^1/2 


^  (  E2  _ 


.,2-2-12/  e2 


,2^1/2 


El{^)  =  ^  (a"  +  +  d")  ,  and  E\{^')  ^  ^  b^  -  d^)  , 

where 

=  [(a^  -  b^f  +  a^bm- 

Products  of  these  form  the  first  few  external  ellipsoidal  harmonics,  which  can  often  be 
scaled  to  give  simple  multiples  of  coordinates  x,  y  and  z; 

X  =  E}(^i)  Em  Em  , 
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J  = 


Eti^,)Eti^2)Eti^2) 


a  b 


(15) 


yz  =  £2(^1) Eldi) £|«3)  , 

ab^ya^  -  b^ 


p2/£  \  172/  £  \  172/  £  \ 

,  ;-o-  -/o  '^2(^1)  £2(^2)  ^2(^3)  2 

-  b^ 


absa^  -  b^ 

However,  external  harmonics  based  on  eI  and  E\  are  weighted  sums  of  x^,y^,  and  ; 


eI{^i)eI{^2)eI{^2)  =  +  a  +  a^z^ +  a,  ,  (16) 


where 


=  (a^  +  b^  +  (f)  {a^  -  2b^  +  (f)l9  , 
ay  =  {a^  +  b^  +  (f){b^-la^  +  (f)l9  ,  (17) 

a,  =  {a^  +  b^  +  (f)  {a^  +  b^-  1(f)  19  , 


a,  =  {f  +  f  +  f)  [  {f  -  ff  -  Ia2b2  +  if  +  f-  f)f]  /  27  , 


and 


E\i^,)E\i^2)Eii^2)  =  a\ 


2  ,  / 

X  +  a 


2  ,  / 

,y  +  a 


2  ,  / 

■z  +  a 


(18) 


where  a\,  a'y,  a\,  and  a'^,  are  given  by  formulas  (17)  with  a'  substituted  for  a  and  -f 
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'y 

substituted  for  d  .  Lame  functions  of  the  second  kind  are  more  complicated; 

FP(^)  =  (2m  +  l)FP(^)  I  (^2_^2)1/2(^2_^2)1/2[£P(^)]2  <^^9) 

=  (2m  +  l)FP(^)  j  (1  _  «2_^2)1/2(1  _  b^x^nEUyx)f  ■ 

The  individual  can,  in  principal,  be  expressed  in  terms  of  elementary  elliptic  inte¬ 

grals.  Some  are  given  in  appendix. 

For  an  external  source  potential  of  a  single  ellipsoidal  harmonic 
E^{^2)  Emi^j,) ,  solution  of  both  limiting  case  boundary  problems  results  in  poten¬ 
tials  within  the  ellipsoid  of  the  same  functional  form  [Equations  (11)],  For  uniform 
source  magnetic  fields  aligned  with  one  of  the  ellipsoid  axes  the  limiting  case  boundary 
problem  solution  has  the  internal  magnetic  field  aligned  with  the  same  axis.  For  the  low 
frequency  limit  boundary  problem,  the  ‘anomalous’  fields  outside  the  ellipsoid  due  to  the 
presence  of  the  ellipsoid  (the  /?o^m  ter^n  in  Equations  (11)),  are  simply  the  magnetic 
fields  of  the  induced  static  magnetization  inside  the  ellipsoid  Mj  =  (jui  -  //o)Vd.  Given 
Equations  (12)  and  (13),  the  anomalous  fields  for  the  high  frequency  limit  boundary  prob¬ 
lem  are  simply  /?olffl=oo/  /^ol®=o  times  these. 

A  Simplified  Higher  Order  Sphere  Ellipsoid  Approximation 

For  a  single  external  ellipsoidal  harmonic  source  field,  the  solution  on  the  inside  of 
an  ellipsoid  has  the  same  functional  form  as  the  source  field,  for  both  low  frequency  and 
high  frequency  limit  solutions,  hence  one  is  tempted  to  make  a  similar  approximation  to 
the  sphere- spheroid  approximation  for  higher  order  ellipsoidal  harmonic  sources. 
Because  of  symmetry,  none  of  the  higher  order  ellipsoidal  harmonics  induce  any  net 
dipole  moment  in  an  ellipsoid,  so  their  effects  first  show  up  in  the  induced  quadrupole 
moments  of  the  ellipsoid.  In  the  case  of  the  three  linearly  polarized  source  fields  (El Ely 
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and  H^),  the  limiting  ellipsoidal  boundary  value  problems  result  in  uniform  magnetic 
dipole  distributions  in  the  direction  of  the  inducing  field  throughout  the  ellipsoid,  with  net 
magnetic  dipole  density  rates  of 

=  m^it)IV  ,  py  =  my{t)IV  ,  p^  =  m^{t)IV  (20) 

per  unit  nominal  inducing  field  in  the  x,  y,  and  z  directions  respectively,  and 
V  =  47rab'c'/3  is  the  volume  of  the  ellipsoid.  The  computationally  simplest  approxima¬ 
tion  for  the  response  to  higher  order  ellipsoidal  harmonic  source  fields,  is  that  the  ellip¬ 
soid  response  is  proportional  to  the  local  magnetic  field  with  the  same  proportionality 
factors  p^,  Py,  and  p^,  as  it  is  for  the  local  magnetic  fields  in  these  directions  for  linearly 
polarized  source  fields.  Under  this  assumption,  because  of  symmetry,  Eli^)  (p  =  yz), 
eI(^)  (P  =  xz),  E\{^)  (P  =  xy),  £'o('?)’  E\{^)  sources  result  in  no  net  induced  dipole 

moment,  and  E\{^)  (p  =  x),  E^(^)  (p  =  y),  and  E^(^)  (p  =  z)  sources  result  in  no  net 
induced  quadrupole  moments.  The  same  assumption  leads  to  induced  quadrupole 
moment  rates  of 


Qy,(t)  =  Hy,  {clffPy  +  blffP,)  y/5  ,  (21) 

2«(0  =  iclffP^  +  alffP^)  y/5  , 

2xy(0  =  H^y  (blffp^  +  alffPy)  y/5  , 

QM  =  2H,AolffPu)VI5  , 

for  source  fields  H  =  (mw  -i-  wu),  with  u  and  w  being  any  of  v,  y,  and  z,  u  and  w 
being  unit  vectors  in  the  corresponding  directions,  and,  according  to  the  choice  of  u,  o^jf 
being  the  corresponding  choice  from  a^jf,  bgjf,  and  Cgjf.  Defining  the  quadrupole 
polarizability  as  being  the  quadrupole  moment  per  for  a  H  =  (mw  -F  wu) 
source  field,  then  =  QuwIHuw  ■  So,  under  the  assumption  of  pointwise  consistency 
between  the  ratios  of  primary  fields  in  the  ellipsoid  axis  direction  and  induced  dipole 
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moment  densities  in  those  directions  within  the  ellipsoid,  one  can  write  the  ellipsoid 
quadrupole  polarizability  rates  in  terms  of  the  ellipsoid  dipole  polarizability  rates,  and  the 
effective  sphere  dimensions  a^jf,  b^jf,  and  Cgjf  for  the  ellipsoid; 

qy^(t)  =  (clffmy(t)  +  blffm,(t))l5  ,  (22) 

<?«(0  =  +  alfftn^it))  15  , 

qxyit)  =  (blffin^it)  +  alffiUyit))  15  , 

=  2og^m„(0 /5  , 

where  u,  w,  and  o^ff  are  as  before.  Quadrupole  polarizability  (rates)  q  have  the  dimen- 
sions  of  Amp-m  /(s-Amp/m  )  =  m  /s.  If  the  equivalent  dipole  polarizabilities  of  an  ellip¬ 
soid  are  known  by  any  means,  quadrupole  polarizabilities  can  be  computed  (for  assumed 
or  known  ellipsoid  dimensions  d,  b'  and  c)  using  Equations  (22),  and  corrections  can  be 
made  for  the  effects  of  quadrupole  moments  due  to  source  field  gradients. 

Under  the  above  assumption,  that  quadrupole  responses  within  the  ellipsoid  are 
proportional  to  the  local  source  magnetic  field  with  the  same  proportionalities  as  the 
dipole  polarizability  responses,  the  resulting  quadrupole  moments  are  linear  in  the  princi¬ 
pal  polarizabilities,  and  the  fields  of  those  moments  are  similarly  linear  in  the  principal 
polarizabilities.  Letting  be  the  anomalous  field  measured  at  the  /’th  sensor  for  the 
j’th  source,  be  the  fields  at  the  i’th  receiver  of  a  dipole  located  at  the  object  cen¬ 
ter,  be  the  fields  at  the  i’th  receiver  of  a  e^e„  quadrupole  located  at  the  object  cen¬ 

ter  (where  ej  =  x,  §2  =  y,  §3  =  z),  be  the  fields  of  the  /th  source  at  signal  maximum, 
at  the  object  center,  be  the  gradient  in  object  principal  coordinates  of  the  fields  of 

the  y’th  source  at  the  object  center,  with  components  and  O  be  the  rotation 

matrix  from  object  principal  coordinates  to  field  coordinates  then  the  anomalous  fields  at 
the  receiver  due  to  the  object  can  be  written  as 
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B 


(a) 


Z^Bf’^^OknPnOlnB 


k,l,n=l,3 


(src) 

Ij 


+  z 

k,l,m,n=\3 


Y^quad)  ^  ^  V7  „  ^2 

U  f^lV  D  J^lj  Pk^l 


/5 


(23) 


where  p^  =  m^{t),  P2  =  my(t),  p^  =  mjp),  =  al^,  =  blj^,  and  r|  =  Expression 
(23)  has  been  simplified  recognizing  that  e„je„  quadrupoles  and  e„e„,  quadrupoles  are 
identical  and  that 


Inversion  for  Ellipsoid  Effective  Semi-axes 

The  anomalous  magnetic  fields  given  by  Equation  (23)  are  linear  in  principal  polar¬ 
izabilities,  granted  that  object  center  location,  orientation  (e.g.,  Euler  angles),  and  effec¬ 
tive  semi- axis  lengths  a^ff,  bgjf,  and  Cgjf  are  known.  We  invert  for  these  parameters  using 
a  general  non-linear  optimization  method  (Smith,  et  al,  1994,  Smith  and  Morrison,  2005), 
at  each  step  fitting  the  principal  polarizabilities  using  linear  methods.  We  minimize  a 
robust  loss  function  of  the  data  residuals  weighted  inversely  by  their  estimated  errors, 
minimizing  squared  weighted  residuals  for  absolute  weighted  residuals  less  than  1  and 
absolute  weighted  residuals  for  absolute  weighted  residuals  greater  than  one  (Huber 
weights,  e.g.,  Huber,  1981). 

The  algorithm  was  tested  on  data  collected  at  the  U.S.  Army  Yuma  Proving  Ground 
calibration  grid,  in  Arizona.  Data  from  vertical  (±90°)  and  dipping  (±45")  81mm  (0.49  m 
long)  and  155mm  (0.87  m  long)  UXO  were  inverted  as  these  ordnance  were  sufficiently 
large  and  shallow  that  source  gradient  effects  were  expected.  Results  using  ellipsoidal 
quadrupole  polarizabilities  are  compared  against  results  from  inversion  for  object  depth 
and  polarizabilities  without  including  ellipsoidal  quadrupole  polarizabilities  in  Tables  I, 
II,  and  III.  As  azimuth  is  undefined  for  a  vertical  object,  azimuth  estimates  for  vertical 
UXO  have  been  omitted. 
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UXO 

Dip  (deg.) 

Depth  (m) 

Est.  Depth  (m) 

Est.  Depth  (m) 

w/o  quadr.  Polariz. 

w/  quadr.  Polariz. 

81mm 

-45 

0.87 

0.80 

0.96 

81mm 

45 

0.87 

0.83 

0.83 

81mm 

-90 

0.95 

0.74 

0.80 

81mm 

90 

0.95 

0.83 

0.92 

155mm 

-45 

1.26 

1.13 

1.14 

155mm 

45 

1.26 

1.09 

1.15 

155mm 

-90 

1.39 

1.22 

1.28 

155mm 

90 

1.39 

1.16 

1.30 

Table  I. 

Estimated  depths  below  the  Berkeley  UXO  Discriminator  (BUD),  for  81mm 

155mm 

UXO  with  and  without  ellipsoid  quadrupole  polarizabilities. 

UXO 

Depth  (m) 

Dip  (deg.) 

Est.  Dip  (deg.) 

Est.  Dip  (deg.) 

w/o  quadr.  Polariz. 

w/  quadr.  Polariz. 

81mm 

0.87 

-45 

-27.2 

-28.0 

81mm 

0.87 

45 

53.3 

52.6 

81mm 

0.95 

-90 

-83.7 

-85.0 

81mm 

0.95 

90 

77.5 

77.2 

155mm 

1.26 

-45 

-22.4 

-57.2 

155mm 

1.26 

45 

23.4 

30.0 

155mm 

1.39 

-90 

-46.9 

-84.0 

155mm 

1.39 

90 

84.1 

85.8 

Table  II.  Estimated  dips  (±180°)  for  81mm  and  155mm  UXO  with  and  without  ellipsoid 
quadrupole  polarizabilities. 
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UXO 

Depth  (m) 

Dip  (deg.) 

Azim.  (deg). 

Est.  Azim.  (deg.) 

w/o  quadr.  Polariz. 

Est.  Azim.  (deg.) 

w/  quadr.  Polariz. 

81mm 

0.87 

-45 

0 

4.3 

6.8 

81mm 

0.87 

45 

0 

-9.3 

-12.9 

155mm 

1.26 

-45 

0 

36.8 

17.8 

155mm 

1.26 

45 

0 

-6.1 

-5.7 

Table  III.  Estimated  azimuths  (±180")  for  81mm  and  155mm  UXO  with  and  without 
ellipsoid  quadrupole  polarizabilities. 


Depths  estimated  without  ellipsoid  quadrupole  polarizabilities  have  a  root  mean  squared 
(rms)  error  of  0.15  m,  those  estimated  with  ellipsoid  quadrupole  polarizabilities  have  0.10 
m  rms  error.  Dip  estimates  have  a  20. 8"  rms  error  without  ellipsoid  quadrupole  polariz¬ 
abilities,  and  a  11.0"  rms  error  with  ellipsoid  quadrupole  polarizabilities.  Azimuth  esti¬ 
mates  have  19. 3"  and  11. 9"  rms  errors  without  and  with  ellipsoid  quadrupole  polarizabili¬ 
ties,  respectively.  In  summary,  using  ellipsoid  polarizabilities  when  estimating  position, 
dip,  azimuth,  and  polarizabilities,  for  objects  where  source  gradients  are  significant  along 
the  length  of  the  object  improves  the  estimates  of  position,  dip  and  azimuth. 

Comparison  of  polarizability  estimates  is  more  difficult  as  true  uniform  field  polar¬ 
izabilities  are  not  known  a  priori,  but  are  estimated  from  measured  data.  Polarizability 
estimates  made  without  ellipsoid  quadrupole  polarizabilities  are  shown  in  Figure  1,  for 
the  8 1  mm  UXO  in  the  four  orientations  of  Tables  I  and  II,  while  polarizability  estimates 
made  with  ellipsoid  quadrupole  polarizabilities  are  shown  in  Figure  2.  Similarly,  polariz¬ 
ability  estimates  made  without  ellipsoid  quadrupole  for  the  155  mm  UXO  in  the  same 
orientations  are  show  in  Figure  3,  those  made  with  ellipsoid  quadrupole  polarizabilities 
are  shown  in  Figure  4.  In  Figures  1  and  2,  and  of  Figures  3  and  4,  the  individual  curves 
are  more  easily  recognized  as  axially  symmetric  objects  in  the  curves  calculated  without 
ellipsoid  quadrupole  polarizabilities  (Figures  1  and  3),  as  there,  the  two  minor 
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polarizabilities  match  each  other  more  closely.  However,  the  polarizability  estimates 
made  with  ellipsoid  quadrupole  polarizabilities  match  each  other  better  between  upward 
and  downward  orientations  (negative  and  positive  dips)  for  both  the  81mm  UXO  and  the 
155mm  UXO,  than  do  the  polarizability  estimates  made  without  ellipsoid  quadrupole 
polarizabilities,  suggesting  that  it  may  be  easier  to  recognize  the  corrected  curves  as 
being  due  to  the  same  UXO,  so  may  be  preferable  for  UXO  model  identification. 

Conclusion 

The  quadrupole  polarizabilities  we  fit  here  were  based  on  assuming  point  wise  con¬ 
sistency  between  the  ratios  of  primary  fields  in  the  ellipsoid  axis  direction  and  induced 
dipole  moment  densities  in  those  directions  within  the  ellipsoid,  between  the  cases  of 
induction  by  uniform  fields  and  induction  by  first  order  gradient  fields.  For  UXO  with 
significant  source  field  gradients  along  their  lengths,  fitting  ellipsoid  quadrupole  polariz¬ 
abilities  under  that  assumption  simultaneous  with  standard  dipole  polarizabilities  reduced 
the  errors  made  in  estimating  object  center  position  and  object  orientation,  and  resulted  in 
polarizability  curves  that  vary  less  with  object  orientation,  at  the  expense  of  poorer  agree¬ 
ment  between  a  particular  response’s  two  minor  (transverse)  polarizabilities.  A  more 
elaborate  ellipsoid  sphere  approximation  for  the  higher  order  terms,  scaling  sphere 
responses  to  yz,  xz,  xy,  and  -  z^  gradient  fields,  in  a  manner  similar  to  Equations  (6), 
to  approximate  ellipsoid  responses  to  yz,  xz,  xy,  Eq{^),  and  E\{^)  gradients,  may  be  more 
accurate  than  the  ellipsoid  quadrupole  polarizabilities  estimated  using  the  assumption  of 
point  wise  consistency  of  induced  dipole  moment  densities  which  were  described  here. 
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Appendix:  Lame  Functions  of  the  Second  Kind  and  Their  Derivative 

Lame  functions  of  the  second  kind  involve  elliptic  integrals.  Here  pertinent  are 
written  in  terms  of  elliptic  integrals  amenable  to  evaluation  using  published  algorithms 
(e.g.,  Bulirsch,  1965a,  1965b,  1969). 
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Derivative  F'^{p)  can  be  expressed  in  terms  of  derivative  E'^{^)  using  Equation  (19)  and 
the  fundamental  theorem  of  calculus,  yielding 
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Figure  (1).  Polarizability  estimates  made  without  ellipsoid  quadrupole  polarizabilities  for 
a  81  mm  UXO  (top  left)  oriented  vertically,  pointed  up,  centered  0.95  below  the  Berkeley 
UXO  Discriminator  (BUD),  (top  right)  oriented  vertically,  pointed  down,  centered  0.95 
below  BUD,  (bottom  left)  dipping  -45°,  centered  0.87  below  BUD,  (bottom  right) 
dipping  45°,  centered  0.87  below  BUD. 
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Figure  (2).  Polarizability  estimates  made  with  ellipsoid  quadrupole  polarizabilities  for  a 
81  mm  UXO  (top  left)  oriented  vertically,  pointed  up,  centered  0.95  m  below  BUD,  (top 
right)  oriented  vertically,  pointed  down,  centered  0.95  m  below  BUD,  (bottom  left) 
dipping  -45°,  centered  0.87  m  below  BUD,  (bottom  right)  dipping  45°,  centered  0.87  m 
below  BUD. 
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Figure  (3).  Polarizability  estimates  made  without  ellipsoid  quadrupole  polarizabilities  for 
a  155  mm  UXO  (top  left)  oriented  vertically,  pointed  up,  centered  1.39  m  below  BUD, 
(top  right)  oriented  vertically,  pointed  down,  centered  1.39  m  below  BUD,  (bottom  left) 
dipping  -45°,  centered  1.26  m  below  BUD,  (bottom  right)  dipping  45°,  centered  1.26  m 
below  BUD. 
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Figure  (4).  Polarizability  estimates  made  with  ellipsoid  quadrupole  polarizabilities  for  a 
155  mm  UXO  (top  left)  oriented  vertically,  pointed  up,  centered  1.39  m  below  BUD,  (top 
right)  oriented  vertically,  pointed  down,  centered  1.39  m  below  BUD,  (bottom  left) 
dipping  -45°,  centered  1.26  m  below  BUD,  (bottom  right)  dipping  45°,  centered  1.26  m 
below  BUD. 


